Quantum Condensation from a Tailored Exciton Population in a Microcavity 
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An experiment is proposed, on the coherent quantum dynamics of a semiconductor microcavity 
containing quantum dots. Modehng the experiment using a generalized Dicke model, we show that 
a tailored excitation pulse can create an energy- dependent population of excitons, which subse- 
quently evolves to a quantum condensate of excitons and photons. The population is created by 
a generalization of adiabatic rapid passage, and then condenses due to a dynamical analog of the 
BCS instability. 

PACS numbers: 71.36.-l-c, 71.35.Lk, 78.67.Hc, 03.75.Kk, 42.55.Sa 
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There is great interest in the possibility of quantum- 
condensed phases of solid-state quasiparticles, such as 
excitons, polaritons, and magnons. Such phases are char- 
acterized by the presence of a quantum state whose pop- 
ulation scales with the size of the system, and hence is 
much larger than one - macroscopic occupation. This 
is seen in recent experiments on Bose-Einstein condensa- 
tion (BEC) of polaritons and polariton lasing. In these 
experiments [H, [3, S 1 a semiconductor microcavity 
is excited at high energies, and a macroscopic population 
of low-energy polaritons emerges following relaxation and 
inelastic scattering Q. These condensates appear spon- 
taneously, from states without macroscopic occupations. 
This differentiates them from the microcavity paramet- 
ric oscillator experiments 0] , where resonant pumping of 
polaritons leads directly to a macroscopic occupation. 

The aim of this paper is to show how microcavities 
could be used to access condensation, even in the ab- 
sence of relaxation or inelastic scattering. We propose 
an experiment on a microcavity containing an ensemble 
of quantum dots, where the exciton decay times are many 
tens or hundreds of picoseconds [§] . We demonstrate that 
this experiment could be faster than these decay times, 
so that energy relaxation and inelastic scattering would 
be negligible. Nonetheless, we shall show that a conden- 
sate develops. In contrast with a laser, this condensate 
is formed from part-matter, part-light quasiparticles. In 
contrast with the microcavity parametric oscillator, it 
develops from a state with no macroscopic occupations, 
in the absence of the pump laser. And whereas relax- 
ation is essential to obtain an equilibrium BEC or po- 
lariton laser, in our approach condensation occurs due to 
an instability of the coherent quantum dynamics. Our 
proposal implements, in a solid-state system, the type 
of dynamical condensation predicted in quenched Fermi 
gases [i, Him, El- 

The first stage in our proposed experiment is the cre- 
ation of a population of excitons in the quantum dots. 
We propose using a chirped laser pulse, which sweeps up 
through part of the inhomogeneously-broadened exciton 
line. As shown by the demonstration of Rabi oscilla- 



tions [l3| and density-matrix tomography [IJ], excitons 
in quantum dots are discrete two-level systems, which 
can therefore be manipulated using laser pulses. The pro- 
posed pulse implements adiabatic rapid passage, which is 
a well-established technique for populating discrete states 
[lit . It extends the technique, by controlling the pump 
spectrum to create an energy-dependent exciton popula- 
tion (Fig. El). 

The second stage occurs after the pump pulse has 
passed. It is the coherent quantum dynamics of the sys- 
tem, starting from the exciton population created by the 
pump. The pump is chosen such that this population 
is similar to a Fermi distribution, with a sharp upper 
step. The system is described by a model similar to that 
which describes pair condensation in atomic gases and 
superconductors. We therefore expect that a population 
with the form of a Fermi distribution could condense, due 
to a dynamical version of the BCS instability [l^ . 

We now turn to the theoretical demonstration of this 
proposal. For simplicity, we suppose that the pump is 
circularly polarized, so we may consider only one of the 
polarization states of the excitons. We model the quan- 
tum dots as a set of two-level systems, each describing 
the presence or absence of an exciton of the pump po- 
larization in a given localized dot state. These localized 
excitons are coupled to the electromagnetic field by the 
dipole interaction. Since the exciton states in the dots 
are spatially separated, we neglect the non-radiative in- 
teractions between the different two-level systems. The 
appropriate Hamiltonian is then the generalized Dicke 
model [lij . 

We label the dot states with an index i, so that Ei is 
the energy of an exciton in the i*^ dot state; this state 
is localized at r^, with dipole-coupling strength gi. The 
area density of dots is n, so that ii gi = g and Ei = E 
the vacuum Rabi splitting at resonance is 2g^/n. The 
state of a dot is specified by the Bloch vector (cr) , where 

= a^—ia^ is the exciton annihilation operator, and 
the inversion. The inversion is related to the occupation 
of the dot rii by rii — {a^ + l)/2, and is — 1(+1) for an 
unoccupied (occupied) dot. 
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Dicke models have previously been used to describe 
polariton condensation in equilibrium [2, 1^, 17 1, and in 
a dissipative open system [isj. Here we are concerned 
with the opposite limit, of timescales short compared 
with the relaxation times. The dynamics therefore obeys 
the Heisenberg equation. Since we are concerned with 
condensation phenomena, involving large photon num- 
bers, we treat the field classically. However, we retain 
the full quantum dynamics of the dots, and hence the 
possibility of an incoherent population of excitons. In 
this approximation, the Heisenberg equation gives 



ill)] 



gPkAE,9)dEdg + Fk, 



(la) 



(lb) 

*i)k,t(£;,5) =.9E (^k',ta+k'.t(i?,g) 

k' (Ic) 

+ P^,_^AE,gH^',t). 

ipk is the normal-mode amplitude for an electromag- 
netic field mode, with in-plane wavevector k and energy 
Wk, and -Fk(^) is the (classical) pump. Pk,t(-B,5) and 
Dk,t{E,g) are collective variables describing the polar- 
izations and inversions of the dots. They are defined as 



Pk(i?,5)<5i?<55 = 



Dk(i?,g)5M.9 = -^^'(a, 
Av ^ — ' 



2\g-ik.r. 
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where the sums run over the states with E ^ E -\- SE 
and g ^ g + Sg. 

Eqs. (|laffTc|) generalize the Maxwell-Bloch equa- 
tions [31, to allow for the distribution of energies and 
dipole-coupling strengths in the dots. Thus the in- 
version and polarization become distribution functions: 
Dq{E, g)5E5g is the inversion due to states with energies 
E ^ E + 5E and couphngs g ^ g + Sg. 

The fields Pk,t,'4'k,t-: and Pk,t have been normalized 
such that their square magnitudes are particle numbers 
per exciton state. A condensate is characterized by a 
macroscopic occupation number, i.e., one which scales 
with the size of the system. Thus an exciton-photon con- 
densate has at least one Pk,t ^ and one V'k.t N^, 
where N is the total number of dots in the active region 
of the sample. In contrast, in a non-condensed state there 
is at most of order one particle per mode, and the P^ t 
and ipk,t are all ^ N~^/'^. 

The approximation leading to ([T]) is the standard semi- 
classical approximation, used to treat condensates in- 
cluding lasers (i^l, superconductors, BCS superfluids 
@i and polariton condensates [21]. It neglects the 



quantum fluctuations of the electromagnetic field, which 
dominate if the photon number is small, i.e., close to 
threshold with a small number of dots [20|, l21|. Since 
we could have N ~ 10^, the semiclassical approximation 
is in general very well-controlled. However, to obtain 
a correct description of the dynamics, we must supple- 
ment ([1]) with noise terms. Without such terms, the non- 
condensed solution remains a steady-state above thresh- 
old, as in all semiclassical treatments of condensation. 
However, it becomes unstable, and hence is not realized. 
We therefore add a perturbation driving the field into 
(flal) , modeling for example spontaneous emission into the 
cavity modes. The form and strength of this perturba- 
tion does not affect our results: we show results with 
Gaussian white noise, but have obtained similar results 
using a delta-like kick. 

We now specialize ([T]) to develop a simulation of our 
proposed experiment. We introduce imaginary parts to 
Wk to allow for the decay of the microcavity photons, 
with timescales of a few picoseconds [l| . The initial con- 
dition is ((t|) = — 1, and there are many exciton states 
distributed over the active area of the sample. Thus the 
sum in ([3]) is strongly peaked near k = 0, and we approx- 
imate the initial conditions as D]s_^t{E,g) — S]g_DQ^t(E, g). 
The dynamics is then that of a continuous medium due to 
motional narrowing, with the short-range spatial struc- 
ture of the exciton states averaged out on the long scales 
of the photons. 

We consider a plane- wave pump, at a high angle where 
the excitons lie outside the stop-bands of the mirrors. 
The field acting on the dots at this wavevector, kp, may 
then be taken as the driving field. Anticipating our anal- 
ysis of the condensation, we retain only one other mode 
of the field, specifically the confined cavity mode with 
k = 0. This reduces ^ to 



i-0o=wo?Ao + ^ J gPo{E,g)dEdg, 

iPo{E,g) = EPoiE,g) - 2gDo{E,g)^o, 
tPp{E,g) = EPp{E,g) - 2gDo{E,g)Fp, 
zDo{E, g)^g {F;Pj,{E, g) + P;{E, g)Fp 
+%Po{E.9) + P^{E,g)i>o). 
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(4b) 
(4c) 

(4d) 



corresponding to an ensemble of two-level systems, inter- 
acting with two modes of the field. Pp is the polarization 
at the pump wavevector, and Fp the driving field. 

We have simulated our proposed experiment by solv- 
ing ([U numerically, with N — 4500 two-level systems. 
Results are shown in Figs. [HO These results focus on a 
model with a single coupling strength g, and a Gaussian 
distribution of exciton energies with variance . The 
pump is a linearly-chirped Gaussian, 
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FIG. 1: Electromagnetic fields and polarizations as functions 
of time, for the simulation described in the text. Dotted: 
pump field \Fp\ (left axis). Dot-dashed: polarization |Pp| at 
pump wavevector, integrated over dot energies and couplings 
(left axis). Solid: cavity field l-i/iol at k = (right axis). 
Dashed: polarization |Po| at k = (right axis). Inset: spec- 
trum |'(/)o('i') 1^ during the shaded region of the main plot (solid 
curve). The pumped population (dashed curve), exciton en- 
ergy distribution (shading), and energy of the k = cavity 
mode (arrow) are shown for comparison. 

where S = J 2g\ip{t)\dt is the usual pulse area per exci- 
ton state. The pulse time r defines dimensionless times 
and energies, and the zero of energy is chosen at the cen- 
ter of the exciton line. The remaining parameters are 
a = 15h/T, ^{luo) -5h/T, vo = -20h/T, a = S/r^, 
— ^(wo) ~ 1.5/t, g — 13/t and S = Svr. These param- 
eters, with T = 3 ps, are reasonable for a microcavity 
containing interfacial quantum dots. Though there is a 
distribution of g, due to the different sizes of the dot 
states, this does not qualitatively change our results. 

Figs. [TH3] are the key results of this paper, demonstrat- 
ing the scenario outlined in our introduction. Referring 
first to Fig.[l] we see that there can indeed be two stages 
to the dynamics. In the first stage, during the pump 
pulse, Pq and "00 are vanishingly small. Pp does become 
finite, reflecting the fact that the pump laser docs in- 
duce some coherent polarization in the cxcitons. This 
coherence is small ~ Pp w 1/VN « 0.02 - but more im- 
portantly both Pp and Fp have disappeared by the end 
of the pumping stage. We therefore argue that the pump 
produces an incoherent population of excitons. Further- 
more, as shown in Fig. [21 this population has a sharp 
upper step, like that of the Fermi function at a low effec- 
tive temperature T ~ 1/r. 

In the second stage, visible in Fig. [1] we see both V'o 
and Pq building up to values of order N'^. This directly 
demonstrates condensation of excitons and photons. 

To understand the condensation, we consider the dy- 
namics of the pumped population. Linearizing ([T]) gives 
normal modes V'k,* = e^^^-*A^, P^ t = e'^'^'* j^p{E, g). 
Their frequencies Ak = AJ^. -I- iX'-^ obey [13] 

u;^-X^^^J ^^dE, (6) 
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FIG. 2: (Color online) Simulated exciton inversion profile im- 
mediately after pumping, showing the population created by 
the Gaussian pump pulse ([5]) (solid line) , and by a superposi- 
tion of such pulses (inset). The dotted curve is an equilibrium 
exciton distribution with fitted temperature h/{4:.2kT); this is 
0.6 K for T — 3 ps. Arrows mark the energy of the k = cav- 
ity mode. 

where i^{E) = J g^Do{E, g)dg is an optical density of the 
excitons. Applying this result to the state immediately 
after the pump pulse, we find an unstable mode at k = 
0. This instability gives the exponential growth of the 
polarization and field, to macroscopic values, visible in 
Fig. m In fact, ^ predicts instabilities for |k| < kc- The 
fastest-growing instability is at k=0, and as this mode 
grows, it suppresses the gain for the others. It will thus be 
dynamically selected, and wc therefore neglected cavity 
modes with k 7^ in the simulations. 

For the parameters used here, the instability predicted 
by ^ corresponds to the BCS instability in a supercon- 
ductor. This can be seen by considering ^ for a single 
coupling strength, close to an instability. The eigenener- 
gies Ak and growth rates A^ then obey 

a;; = ng\gni\'^)Do{X'^) - 7. (7b) 

(|7b)) describes the growth or decay of the normal mode, 
with the first term the gain/loss from the excitons, and 
the second the loss due to the cavity decay 7 = — S(wk)- 
d7a|) is the Cooper equation of the BCS model [2^. The 
term corresponding to the usual pairing interaction is 
5^/(Ak — i^k), which we recognize as the effective inter- 
action between excitons, mediated by the cavity modes. 
Here it is an attractive interaction, as required for BCS, 
since the excitons lie below the photons. The term corre- 
sponding to the Fermi distribution is the exciton popula- 
tion created by the pump. In a superconductor, there is a 
solution to the Cooper equation below the Fermi energy, 
due to the step in the Fermi distribution; in the same 
way, (|7a| has a solution below the step in the exciton 
occupation. (|7b[) shows that this mode experiences gain, 
and hence can become unstable. 
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FIG. 3: (Color online) Simulated exciton inversion profile dur- 
ing the condensation shown in Fig. [T] 



To confirm the origin of the condensation, we show in 
the inset to Fig. [1] the spectrum of the k = field, and 
in Fig. [3] the evolution of the population during conden- 
sation. As expected, the condensate is at a frequency 
below the step in the exciton population (Fig. [T]); this 
leads to hole-burning there (Fig. [3]). Reducing the center 
frequency of the chirp, j/q, we find that the condensation 
stage in Fig. [T] disappears, as the effective pairing inter- 
action decreases below its critical value. The long-time 
limit of the simulations is then just the population of 
excitons. We have also confirmed that the phase of the 
pump ^ is irrelevant, by plotting the phase associated 
with the power spectrum in the inset to Fig. [T] This plot 
is identical in simulations with different pump phases, 
but the same random noise. 

Since our condensate occurs on timescales short com- 
pared with the relaxation times, it will not be in equilib- 
rium. This leads to phenomena absent for an equilibrium 
condensate. In Fig.[l] for example, we see ringing oscilla- 
tions (corresponding to those predicted for atomic gases 
U , TM y unrelated oscillations have been predicted 
in coherently-driven microcavities (23j). and a slow decay. 
These phenomena will be discussed in a future publica- 
tion. 

It is perhaps surprising that our pump configuration 
could generate an incoherent population, without macro- 
scopic occupation. To explain the mechanism, we note 
first that the driving field dominates during the pump 
stage (Fig. [1]). Thus the pumping can be understood 
in terms of the well-known dynamics of non-interacting 
two- level systems, driven by a chirped laser pulse [l5| . 
After eliminating the time-dependent frequency uj{t) of 
the pump with a unitary transformation, this dynam- 
ics is a precession of the Bloch vectors (cr^) around axes 
Bi — {2giR{t),0, uj{t) — Ei), where R{t) is the pump field 
at the dot. 

Consider first a dot with energy inside the chirp. Pro- 
vided the pump is strong enough, the Bloch vector of 
such a dot adiabatically follows from — k to k, and 
it becomes populated. Furthermore, a dot with energy 



outside the range of the chirp will not respond unless the 
pump is very strong, and hence such a dot will remain un- 
populated. Thus we see that the only dots which could be 
polarized by the pump are those at the edge of the chirp. 
However, there are very few such dots. Furthermore, 
any polarizations induced in dots with different energies 
would have different phases, due to the chirping. Thus, 
as shown in Fig. [U the chirped pump does not induce 
a collective polarization. Noting additionally that there 
are no cavity modes resonant with the pump, we see why 
the pump generates a state with Pk,t,'!/'k,t ^ N~^/'^. 

A fully quantitative description of our proposed ex- 
periment is likely to require a more detailed model of the 
quantum dot states. It would also be desirable to extend 
our model to (a)incorporate relaxation processes, which 
on longer timescales will drive the off-equilibrium con- 
densate towards an equilibrium one and (b)develop 
a full quantum theory, as has been done for photon lasers 
2C I , polariton lasers US] j and equilibrium condensates 
211 . Finally, it would be interesting to compare this 
work with Ref. [2^, which appears to show a type of 
dynamical condensation, under very different conditions. 

Though the parameters used here have been chosen to 
suit interfacial quantum dots 0, [2^ , many other choices 
are possible. The basic requirements are that the pump- 
ing is slow enough to create a controlled inversion pro- 
file within the inhomogeneous line, such that the cou- 
plings can be increased so that the instability occurs af- 
ter the pump pulse, but before the excitons decay. Other 
systems described by generalized Dicke models, such as 
Fermi gases, SK dots in microcavities, or Josephson junc- 
tion arrays, could be considered. 

To conclude, we have proposed and analyzed a new 
approach to quantum condensation in a solid-state sys- 
tem. The key is the chirped pumping ([S]), which we have 
shown can create an energy-dependent population in the 
exciton line. Such a population can condense, even in 
the absence of relaxation or inelastic scattering. Since 
our approach uses the spectrum of the pump to tailor 
the population, it would be possible to pump other ini- 
tial states (inset to Fig. [2]). Thus our technique could be 
used more generally, to explore the quantum dynamics of 
a many-particle system from controlled initial conditions. 
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EP/F040075/1 and EP/C546814/01. 



[1] H. Deng, D. Press, S. Gotzinger, G. S. Solomon, R. Hey, 
K. H. Ploog, and Y. Yamamoto, Phys. Rev. Lett. 97, 
146402 (2006). 

[2] J. Kasprzak, M. Richard, S. Kundermann, A. Baas, 
P. Jeambrun, J. Keeling, F. Marchetti, M. Szymanska, 
R. Andre, J. Staehh, et al.. Nature 443, 409 (2006). 

[3] R. Bahh, V. Hartwell, D. Snoke, L. Pfeiffer, and K. West, 
Science 316, 1007 (2007). 



5 



[4] D. Bajoni, P. Senellart, E. Wertz, I. Sagnes, A. Miard, 
A. Lomaitre, and J. Bloch, Phys. Rev. Lett. 100, 047401 
(2008). 

[5] A. P. D. Love, D. N. Krizhanovskii, D. M. Whittaker, 

R. Bouchekioua, D. Samvitto, S. Al Rizeiqi, M. S. Skol- 

nick, P. R. Eastham, R. Andre, and L. S. Dang, Phys. 

Rev. Lett, (in press). 
[6] T. D. Doan, H. T. Cao, D. B. Tran Thoai, and H. Haug, 

Phys. Rev. B 74, 115316 (2006). 
[7] R. M. Stevenson, V. N. Astratov, M. S. Skolnick, D. M. 

Whittaker, M. Emam-Ismail, A. L Tartakovskii, P. G. 

Savvidis, J. J. Baumberg, and J. S. Roberts, Phys. Rev. 

Lett. 85, 3680 (2000). 
[8] W. Langbein and B. Patton, Phys. Rev. Lett. 95, 017403 

(2005). 

[9] A. V. Andreev, V. Gurarie, and L. Radzihovsky, Phys. 
Rev. Lett. 93, 130402 (2004). 
[10] R. A. Barankov and L. S. Levitov, Phys. Rev. Lett. 93, 
130403 (2004). 

[11] E. A. Yuzbashyan, O. Tsyplyatyev, and B. L. Altshuler, 

Phys. Rev. Lett. 96, 097005 (2006). 
[12] P. R. Eastham, J. Phys.: Condens. Matt. 19, 295210 

(2007). 

[13] A. J. Ramsay, R. S. Kolodka, F. Bello, P. W. Fry, W. K. 
Ng, A. Tahraoui, H. Y. Liu, M. Hopkinson, D. M. Whit- 
taker, A. M. Fox, ct al, Phys. Rev. B 75, 113302 (2007). 

[14] Y. Wu, X. Li, L. M. Duan, D. G. Steel, and D. Gammon, 
Phys. Rev. Lett. 96, 087402 (2006). 



[15] V. S. Mahnovsky and J. L.Krause, Eur. Phys. J. D 14, 

147 (2001). 

[16] J. Keeling, P. R. Eastham, M. H. Szymanska, and P. B. 

Littlewood, Phys. Rev. B 72, 115320 (2005). 
[17] P. R. Eastham and P. B. Littlewood, Phys. Rev. B 64, 

235101 (2001). 

[18] M. H. Szymanska, J. Keeling, and P. B. Littlewood, Phys. 

Rev. B 75, 195331 (2007). 
[19] K. Staliiinas and V. J. Sanchez- Morcillo, Transverse 

Patterns in Nonlinear Optical Resonators, no. 183 in 

Springer Tracts in Modern Physics (Springer, 2003). 
[20] H. Risken (North Holland, Amsterdam, 1970), vol. 8 of 

Progress in Optics, chap. 5, pp. 241 - 294. 
[21] P. R. Eastham and P. B. Littlewood, Phys. Rev. B 73, 

085306 (2006). 

[22] P. G. de Gennes, Superconductivity of Metals and Alloys 

(Westview, 1999). 
[23] L A. Shelykh, M. M. Glazov, D. D. Solnyshkov, N. G. 

Galkin, A. V. Kavokin, and G. Malpuech, Phys. Stat. 

Sol. (c) 2, 768 (2002). 
[24] F. P. Laussy, G. Malpuech, A. Kavokin, and P. Bigen- 

wald, Phys. Rev. Lett. 93, 016402 (2004). 
[25] H. Vinck-Posada, B. A. Rodriguez, P. S. S. Guimaraes, 

A. Cabo, and A. Gonzalez, Phys. Rev. Lett. 98, 167405 

(2007). 

[26] A. G. Steffan and R. T. Phillips, Phys. Stat. Sol. (a) 190, 
541 (2002). 



